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ABSTRACl’ 

New integral formulas to determine geopotential coefficients trom terrestrial gravity and satellite alti- 
metry data are given. The formulas arc based on the integration of data over the non-spherical surface of the 
Earth. The effect of the topography to low degrees and orders of coefficients is estimated numerically, for- 
mulas for the solution of the gravimetric boundary value problem are derived. 


1. Introduction 

The long- wavelength features of the Earth’s gravity field is best determined from dynamic satellite ob- 
servations. For short-wavelength information terrestrial gravity anomalies (Ag) and satellite altimetry data 
play essential roles. A well-known method to determine gcopotential coefficients from Ag is the integral met 
hod used by Rapp e. g. (1977). We present this method in section 2 followed by a new integral method in sec- 
tion 3. 


2. Rapp s Integral Formula 

The determination of the harmonic coefficients according to Rapp requires that the data (Ag) is distri- 
buted continuously all over the surface of the Earth. Also it is assumed that the Earth’s ( disturbing) potentia 
(T) can be developed into a harmonic series convergent down to the Earth’s surface: 

00 Tq . i n /o i \ 

T = r 0 -y X (-) 2 A nm Y nm (<J>,\), (2<1 > 

n = 2 r m = -n 

where 7 is the mean surface gravity of the Earth, (r,4>,\) are the geocentric, spherical coordinates of the 
computation point, Y n m(<t>,X) is a fully normalized spherical harmonic and, for r = r 0 , 

Amn = // T Y mn da, ( 22 > 

4Trr 0 7 a 

where a is the unit sphere. Considering the "boundary condition" in spherical approximation (Heiskanen and 
Moritz 1967, p 87) one obtains the corresponding expansion for the gravity anomaly (Ag) : 


Ag = y 2 (n ~ 1)( — ) n + 2 5) Anm Ynm (<)>A) j 
n = 2 r m = -n 

from which formula the coefficients are given for r = r 0 : 


(2.3) 


Anm = ff Ag Y nm da . ( 2A) 

4Ti7(n-l) a 

In principle r 0 is arbitrary, but is usually chosen as the equatorial radius of the selected reference ellip- 
soid. A practical obstacle in applying the integral formulas (2.2) and (2.4) is the necessity of reducing T and 
Ag to the sphere of radius r Q . For (2.4) Rapp (1984) used the Taylor series 


3Ag 1 d 2 Ag 2 

Agrcduccd — Agobscrvcd Ah (Ah) - ... , 

ah 2 ah 2 


(2.5) 


where Ah = r - r () , and the derivatives were estimated from available scries of harmonic coelficienls. 
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3. The New Integral approach 


Rewriting Tormula (2.1) on the form 

T = R 7 1 ( — ) n + ’ 1 CnmY n m(4>A), 

n = 2 r m = - n 

where R is the radius of die minimum bounding sphere of the Earth, and considering Green s second identity 
for the functions T and U (Hciskanen and Moritz, 1967, p. 11): 

ST 6U . ~ 

f/J (UAT-TAU)dv = // ( U T -) dS, 

v S Sn’ 6n’ 

where v is the volume bounded by the surface S, which consists of the surface of the Earth (E) and the sur- 
face of the bounding sphere of radius R, and n’ is the external normal to S with respect to v, and setting U to 


U = { ( — ) n+ 1 + an ( — ) n } Ynm(4>A)> 

r R 


(3.3) 


where a n is a non-zero constant with respect to position, the following general integral formula is derived in 
Sjoberg (1988): 

1 — 1 ~ 8 , 2 


Cnm — 


4TTR~7(2n 4- l)an 




[U { Ag — T } + T ] r da, 

y Sr Sr 


(3.4) 


where Ag = Ag - y(£ tan + t\ tan 32 ), pi and p2 are the terrain inclinations, and t, and p are the corre- 
sponding deflections of the vertical. The integral is taken over the Earth’s surface (E). Choosing r-r 0 + Ah 
and a n = (R/r 0 ) 2 " + 1 (n - l)/(n + 2) one obtains (to order Ah/r 0 ): 


ELEVATION EFFECT Case I 
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R 

Anm = ( ) Cnm Anm "hAAnm, 

To 


where A n m° is given by (2.4) and 

n + 2 

AAnm = // T All Vnni d(J\ 

4TryRr 0 K 


(3-5) 


(3.6) 


The percentage terrain corrections by degree determined by formula (3.6) are illustrated in Hg. 1. Deta- 
ils on the computations are given in Sjoberg and Fan (1988). For comparison the Figure includes also terrain 
corrections from Rapp (1984). 


Choosing a n = - (R/r 0 ) 2n + 1 one arrives at formula (3.5), but now with A n m° given by (2.2) and 
1 R _ 2T Ah 

AAnm = (“) // (Ag + *“ ) — Y n m(4>A) dfT - 

4 Try r 0 cr R r 0 

Numerical computations show that this correction is within 2% for n < 50. Again we refer to Sjoberg 
and Fan (1988) for details. 


4. Solutions to Boundary Value Problems 

Inserting (3.4) into (3.1) and assuming that summation and integration may change order one arrives at 
the following solutions for the height anomaly £p = Tp/7 : 


Case I: a n = (R/r 0 ) 2n + 1 (n-l)/(n + 2) : 

R _ 

£p = JJ ( { S(i|ipQ,rQ,rp) - A(^PQ,rQ,rp) }AgQ + K(i|/pQ,rQ,rp) Tq ] do-Q, (4.1) 

4tt7 o 

where 

tq oo 2n + 1 tq 

S(i|>PQ,rQ,rp) = — £ (— ) p n (cos 4>PQ) 

R n = 2 n — 1 rp 

tq oo n + 2 ro - rQ 

A(^PQTOTP) = — 2 { 1 - (— ) 1 } ( — ) p n (cos 4 /pq) 

R n — 2 n — 1 tq rp 

and 

K(v|ipQ,rQ,rp) = — S (n +2) (— ) n 1 { 1 - (— ) } P n (cos i|ipq) 

R n = 2 rp tq 

Case II: a n = -(R/r G ) 2n + 1 : 

1 _ 

O' = ff [ { P(>W , Qi r Q,rp)-L(<)TO,rQ,ri>) } Tq + M(»]/pQ,rQ,rp) tq Ago ] dciQ , (4.2) 

4Try <j 
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where 

00 r 0 

P(>l'PO,ro,rp) = X(2n+1)(— ) n + 1 Pn(cos^>o) 
n =2 rp 


L(ipPO,rQ,rp) 


= X(n-l) A n + 1 {l- (-) 2n + 1 }P n (cos^PO) 
n = 2 rp tq 


and 


M(i|»pQ,rQ,rp) 


_ £ ( )" + 1 { i . (— ) 2n + 1 } P n (cos ippo). 
n = 2 rp r Q 


l„ ,bc limil ,o - ro Che formulas (4.1) and (4.2) approach (.he emended) Sloke’s formula and the Bruns 
Poisson’s formula, respectively. 


5. Discussion 

As Figure 1 shows the terrain corrections of the new method (3.6) are more significant with mere^ing 

S3 -^SSSS^SSSaSaSa 

from surface data. . 

no. G-GU 4071-302. 
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